A quasi-Hermitian pseudopotential for higher partial wave scattering 
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We formulate a new quasi-Hermitian delta-shell pseudopotential for higher partial wave scatter- 
ing, and show that any such potential must have an energy-dependent regularization. The quasi- 
Hermiticity of the Hamiltonian leads to a complete set of biorthogonal wave functions that can be 
used as a basis to expand and diagonalize other two-body Hamiltonians. We demonstrate this pro- 
cedure for the case of ultracold atoms in a polarization-gradient optical lattice, interacting pairwise 
when two atoms are transported together from separated lattice sites. Here the pseudopotential 
depends explicitly on the trapping potential. Additionally, we calculate the location of trap-induced 
resonances for higher partial waves, which occur when a molecular eigenstate is shifted to resonance 
with a trap eigenstate. We verify the accuracy of the pseudopotential approach using a toy model 
in which a square well acts as the true interaction potential, and see excellent agreement. 

PACS numbers: 34.20.Cf, 34.50.-s, 32.80.Pj 



I. INTRODUCTION 

The behavior of trapped ultracold neutral atoms has 
attracted considerable interest in both quantum con- 
densed matter physics and quantum information science. 
In these disciplines, ultracold atoms in tightly confining 
traps are an attractive platform given the richness of the 
systems and the ability to coherently control a wide va- 
riety of parameters in a precise and well understood way. 
From a condensed matter physics perspective, ultracold 
atoms trapped in optical lattices provide an ideal sys- 
tem for studying tight-binding models such as the Hub- 
bard model [1-3]. In the context of quantum informa- 
tion, optical trap-arrays provide a scalable platform for 
storing many qubits, with parallel operations, applica- 
ble to quantum-cellular automata [4, 5] , or more general 
quantum circuit operations [6], and fault-tolerance via 
topological encoding [7]. 

At the fundamental level, quantum many-body sys- 
tems are governed by their two-body interactions. In the 
case of ultracold neutral atoms, these consist of collisions 
determined by the diatomic molecular interaction poten- 
tial. Coherences between unbound scattering atoms and 
bound molecules provide new handles for tuning interac- 
tions and exploring heretofore uncontrollable many-body 
phenomenen such as superfluid to Mott insulator transi- 
tions [8] and the BCS to BEC crossover [9]. The pro- 
duction of ultracold molecules also opens the door to 
new studies in quantum chemistry where nuclear mo- 
tion is fundamentally quantum and controllable. Addi- 
tional novel phenomena arise in tightly confining traps 
due to reduced dimensionality including the Tonks Gas 
for strongly interacting ID Bose systems [10] and trap- 
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induced resonances [11, 12]. 

For ultracold collisions, the Wigner-threshold law dic- 
tates that the lowest nonvanishing partial wave domi- 
nates the interaction. Unless ruled out by exchange sym- 
metry, s-wave collisions dominate and are generally char- 
acterized by one parameter, the scattering length a. This 
can be codified into a contact "pseudopotential" , whose 
scattering phase shift in the Wigner-threshold regime 
equals that of the true molecular potential. The use 
of the Fermi pseudopotential well-describes the proper- 
ties of degenerate quantum gases, such as the mean-field 
effects predicted by the Gross-Pitaevskii equation [13]. 
The case of two atoms in a trap can be treated with great 
efficacy by the same pseudopotential as for atoms scat- 
tering in free space. The contact potential effectively sets 
the boundary condition for the relative coordinate wave 
function at short range [13-15] whereas the trap sets the 
long-range boundary conditions. The resulting two-atom 
problem determines a bound-state energy spectrum and 
corresponding set of eigenfunctions that can serve as a 
basis for the studies of many-body systems with pairwise 
interactions in a trap [11-13, 15, 16]. 

Although s-wave collisions typically dominate at ul- 
tracold temperatures, in some cases higher partial waves 
play an important role in the interaction. For exam- 
ple, in the case of identical fermions, s-wave scattering 
is forbidden and p-wave scattering dominates. For the 
boson 133 Cs there are anomalous collision properties due 
to bound state of the dimer near zero energy leading to 
large spin relaxation through coupling to higher partial 
waves [17]. The generalization to higher partial wave 
scattering has recently been revisited by numerous au- 
thors. A long-standing error in the Huang pseudopoten- 
tial for arbitrary partial waves [18] has been corrected 
[19-23]. Our solution [19], based on a delta-shell poten- 
tial at finite range, allows for arbitrary multipole distri- 
butions associated with higher-order partial waves. This 
has been used to calculate the level shifts of the eigen- 
states of trapped atoms in the context of optical lattices 
[11, 24]. 
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For s-waves, the pseudopotential also gives rise to nor- 
malizable wave functions that can be used to diagonal- 
ize other Hamiltonians describing interactions between 
trapped atoms, e.g. two atoms in separated traps of a 
polarization-gradient lattice [11, 24]. Such a system ex- 
hibits new trap-induced resonances, which occur when 
a molecular bound state of the interaction potential is 
shifted to the same energy as a vibrational state of the 
trap [11]. These trap-induced resonances offer a new 
handle for the control of interactions between ultracold 
atoms, and for the design of quantum logic gates. In 
other examples, generalized pseudopotentials have been 
used to study higher-order partial wave interactions in 
anisotropic traps, including the approach to quasi-2D 
gases [21, 25]. 

In this article, we introduce a new quasi-Hcrmitian 
higher-partial wave pseudopotential that allows us to 
construct a complete biorthonormal set of cigenfunctions. 
This is essential for the application of the pseudopo- 
tential to trapped atoms, and in particular to the cor- 
rect description of higher partial wave trap-induced res- 
onances. Whereas our previously derived pseudopoten- 
tial [19] based on a delta shell allows a correct derivation 
of the eigenspectrum of interacting atoms, as reviewed in 
Sec. II, the eigenfunctions for this non-Hermitian pseu- 
dopotential do not form an orthonormal set. We show 
that even in the Wigner regime with constant higher 
partial wave scattering lengths, we need an energy- 
dependent regularization to correct this problem. In Sec. 
Ill we construct the appropriate quasi-Hcrmitian poten- 
tial. The time-independent Schrodingcr equation is then 
solved self-consistently, giving rise to a biorthonormal 
set of wave functions, as discussed in Sec. IV. We ap- 
ply our solution to the study of trap-induced resonances 
in polarization-gradient optical lattices, generalized to 
higher-partial wave scattering in Sec.V. Our results are 
tested using a toy model, where a square well acts as the 
true interaction potential, and find excellent agreement. 
Section VI summarizes our results. 



II. REVIEW OF THE PSEUDOPOTENTIAL 
CONSTRUCTION 

Ultracold atoms can have a very long dcBroglie wave 
length, on the order of microns, whereas the van der 
Waals interaction for alkali atoms is on the order of 1-10 
A[26]. For scattering from a short-range, radially sym- 
metric potential, the effect of the interaction is to phase 
shift the radial wave function so that it takes the asymp- 
totic form 

Fi>k(r) = A t (j,(fcr) - tan(J,(fc))n,(fcr)) , (i) 

where 5i is the scattering phase shift, fc is the wave num- 
ber associated with the relative coordinate kinetic energy, 
and ji, m, are the spherical Bessel and Neumann func- 
tions, respectively. It is possible to replace the true in- 
teraction potential with a simpler potential, as long as it 



gives the same asymptotic wave function. For s-waves in 
the Wigner threshold regime, where the scattering phase 
shift is related to the scattering length by tan6 w — ka, 
this can be done using a regularized Fermi's pseudopo- 
tential, defined as 



V (r) = 2™«5( 3 )(r)-r. (2) 



Here and throughout, we take units with h = fi = 1, 
where fi is the reduced mass of the two-atom system. 
The first order derivative serves to regularize the pseu- 
dopotential by removing the singular contributions at 
short range, allowing it to act on arbitrary wave func- 
tions, including those that diverge at the origin, e.g., 
the Neumann functions n; in Eq. (1). The only param- 
eter needed for reproducing the actual asymptotic wave 
function is the zero-energy scattering length, a, given by 
a = - linu^o [tan {S )/k]. 

For two atoms in a tightly confining trap, the interac- 
tion causes a shift in the unperturbed vibrational states. 
There are three length scales that govern the physics: the 
range of the interaction potential, (3§ = (mCs/H 2 ) 1 / 4 for 
van der Waals interactions, the scattering length, a, and 
the width of the unperturbed ground vibrational state 
in the trap, tq. These length scales are often very dif- 
ferent, making it difficult to find a numerical solution 
to the Schrodinger equation for geometries that are not 
separable into ID solutions. A solution can be found 
by modeling the net effect of the interaction by the same 
pseudopotential as in free space if the confining trap does 
not significantly distort the short-range interaction po- 
tential, true when f3§ <C ro, which is usually the case. 
This approach has the advantage of separating the short 
range behavior from the long range solutions. 

In a seminal paper, Busch et al. [14] solved for the 
new energy eigenfunctions and values for two atoms in 
a harmonic trap interacting via s-wave collisions, mod- 
eled using the regularized Fermi pseudopotential, pa- 
rameterized by the free-space scattering length a. The 
solution is found in the transcendental equation, a = 
T(-v - l/2)/2r(-z/) with E = 2v + 3/2, having chosen 
units with the trap fre quency u = 1 so that all lengths 
are measured in tq = \Jh/ (jlj and energy in units of fuv. 
As long as the vibrational energy is within the Wigner 
threshold regime, guaranteed when a <C r , this accu- 
rately captures the perturbed spectrum. Out of that 
regime, one employs a free-space pseudopotential, pa- 
rameterized by an energy- dependent scattering length, 
a — > a(E) = — tan<5o(fc)/fc, where E — k 2 /2. Since now 
the strength of the pseudopotential depends on the vi- 
brational energy, whose value we seek, the Schrodingcr 
equation must be solved self consistently [13, 15]. One 
solves for the energy as a function of scattering length 
E(a) (the "Busch solutions") and separately a(E) in a 
free space scattering calculation. This leads to the fol- 



3 



lowing set of equations, 

tan (So (A:)) _ lT{-v - 1/2) 
fc ~ 2 r(-i/) : 

i/c 2 = 2^ + 3/2. 



(3a) 
(3b) 



These solutions show excellent agreement with full scat- 
tering calculations [15]. 

For higher partial waves, this construction is more sub- 
tle. The potential is separable into a radial and an an- 
gular part, with 



V(r)= Vl (r)Y lm (f), 



(4) 



where Y/ m are the spherical harmonics. A contact poten- 
tial at the origin cannot support higher-order multipoles, 
I > 0. Thus, no pscudopotential for higher partial waves 
can be constructed exactly at the origin. Derevianko [27] 
solved this by employing the appropriate Green's identity 
in order to relate a surface integral at a finite radius to 
the pseudopotcntial. Our solution explicitly defines the 
pscudopotential as the limiting form of a finite radius 5- 
shell [19]. Defining higher-order scattering lengths by the 
Wigner threshold law, tant^(fc) w — (fca;) 2 ' +1 , the <$-shell 
pseudopotcntial takes the form, 



Vi{r) = lim a t 



21+ 



l 6(r - r s ) 



rl +2 



D, 



where 



D, = 



(21 + 1)!! 1 d 2l+1 
2(2/)!! ^i+2dr 2l + 1 '' 



(5) 



(6) 



is the regularization operator. For s-waves, one recovers 
the familiar regularized Fermi pseudopotential. 

As demonstrated in [19], the generalization of the 
Busch solutions to arbitrary partial-wave scattering 
follows from Eq. (5), parameterized by the energy- 
dependent /-wave scattering phase shift in free space ac- 
cording to the substitution a 2l+1 — ► — tan<5;(fc)/fc 2 ' +1 = 
Pi(k) (see also [28] for a quantum-defect theory ap- 
proach). An ansatz for the radial wave function of 

the form F(r) — r L e~ r / 2 w(r) transforms the piecewise 
quadratic radial equation into Kummer's equation, with 
confluent hypergeometric functions M and U as solutions 
[24, 29]. Inside the shell, the solution must be the Rum- 
mer function that is regular at the origin, while outside, 
the boundary equation at infinity requires the solution 
which asymptotically approaches zero, 

w <(r) = Bi M(-v,l + 3/2,r 2 ) r<r s , (7a) 
w > (r) = Ai U(—v, / + 3/2, r 2 ) r > r Sl (7b) 

where the energy eigenvalue is£' = 2^ + / + 3/2in har- 
monic oscillator units. Note, when v is an integer, the so- 
lution everywhere is the familiar Laguerre polynomial for 
an unperturbed isotropic harmonic oscillator; nonintegcr 
v arise from the pseudopential, which shifts the energies 



of the trap eigenfunctions. The boundary conditions for 
the wave function across the delta-shell, 



1 UF< 

2 \~dr~ 



F<(r s )=F>(r s ), 
dF>' 



dr 



= af +1 D l F>{r s ) 1 



(8) 



determine the wave functions and a transcendental equa- 
tion for the eigenvalue v, 



tan 5i(k) 
k 2l+1 



(-1)' 



r(/ + 3/2) 



n 2 



(2/ + 1)!! 



r(-i/-/-i/2) 



(9) 



The generalized, self-consistent Busch eigenvalues are 
those for which E = 2v + I + 3/2 = fc 2 /2. 



III. A QUASI-HERMITIAN 
PSEUDOPOTENTIAL 

Whereas the higher partial-wave pseudopotential sug- 
gested in [19] allowed us to derive the generalized eigen- 
value spectrum for trapped atoms, it fails, as we will 
show in the following, in the construction of a complete 
set of orthonormal eigenfunctions. This task is stymied 
by the regularization operator of the pseudopotential Eq. 
(6) which contains odd powers of derivative operators, 
rendering it non-Hermitian for all but s- waves [32]. Its 
eigenfunctions do not form an orthonormal set, as can be 
seen by considering the limit of the eigenfunctions when 
the shell radius r s of the potential is taken to zero. Ac- 
cording to Eq. (5), the pscudopotential for p- waves has 
the form 



lim c4r<5(r • 



,d 3 2 



(10) 



where C is a constant. For / > 0, the reduced radial 
wave functions outside the shell, rF > (r) in Eq. (7), di- 
verge as r^', for r s — > 0. Therefore, most of the proba- 
bility is centered around r = r s . Normalizing these wave 
functions makes the oscillatory tail negligible as r s — > 0, 
and hence, the integral over wave functions associated 
with different energies cannot be zero because the pro- 
nounced peak common to all wave functions cannot be 
canceled. In contrast, for s-waves the outside reduced ra- 
dial functions rF y (r) go to a finite value at r — ► 0, and 
hence, we can take the limit of zero shell radius with- 
out obtaining this pronounced peak. In fact, in the limit 
r s — > the pseudopotential for s-waves has an orthonor- 
mal, complete set of eigenvectors with real eigenvalues 
and is therefore Hcrmitian. For a more rigorous proof, 
see Appendix. 

We can obtain a complete set of eigenfunctions asso- 
ciated with a pseudopotential for higher partial waves in 
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the following way. Since the Hamiltonian has real eigen- 
values, it is possible to construct a quasi- Hermitian po- 
tential, with a biorthonormal set of eigenvectors, 



if = ]T I, 



(ii) 



such that H\V n ) = E n \$ n ), ($ n \Hi = (® n \E n , and 
(^ n \^m) = &nm [30]. This, however, requires a new reg- 
ularization since a delta-potential with an odd number of 
derivatives beyond first order does not possess a quasi- 
Hcrmitian adjoint with finite extension. Projecting onto 
the partial wave I, we take a new ansatz for the radial 
dependence of the pseudopotential restricted to a regu- 
larization with at most a first order derivative, 



Vi(r) = ci tanSi(k)S(r - r s ) 



d 
or 



(12) 



where C\ and c 2 are constants to be determined and Si(k) 
is the scattering phase shift at kinetic energy Ek = k 2 /2. 

We begin by constructing a higher partial wave pseu- 
dopotential for atoms scattering in free space, where the 
radial wave function corresponding to the eigenstate of 
the pseudopotential takes the form Eq. (1) for r > r s and 
-^i < ( r ) = Biji(kr) for r < r s . The constant c 2 should be 
chosen in such a way as to remove the irregular compo- 
nent of the wave function outside the shell, 



d 

C2 + 7T 
or 



ni{r) = 0. 



(13) 



Additionally, we must satisfy the boundary condition 
given by integration of the Schrodinger equation over the 
shell radius, 



cic 2 tan(<5 / (fc))F > 
d 



+ 



dtan^fc)) j^F> 



0. (14) 



These conditions can be achieved only with a regulariza- 
tion that depends explicitly on the kinetic energy of the 
colliding particles. As we have already considered the ex- 
tension of the pseudopotential with an energy-dependent 
scattering length to account for situations beyond the 
Wigncr-thrcshold regime, the extension to an energy- 
dependent regularization is natural. Though not unique, 
one possible choice for c\ and c 2 is 



ci 



ni{kr s ) 



C2 = - 



n l( kr s) 

ni{kr s )' 



(15) 



For small shell radii, and for s-wave scattering (1 = 0), 
c\ — > —1/kr, c 2 — ► 1/r. The resulting pseudopoten- 
tial is seen to recover the familiar Fermi pseudopoten- 
tial for r s — > in the Wigner threshold regime with 
tan(<5 (fc)) ->■ ~ka with 5(r - r s ) -> ^r 2 5^ (r). Note 



that given the wave functions' dependence on the kinetic 
energy, the pseudopotential for higher partial waves can- 
not be reduced to an energy-independent form, even in 
the Wigner-threshold regime. As a result we must solve 
the eigenvalue equation self-consistently, as will be dis- 
cussed in Sec. IV. 

This quasi-Hermitian pseudopotential can be extended 
to the case of interacting atoms in traps. The radial wave 
functions take the canonical form that generalizes Eq. 
(1), 



F£ E {r) = Bif,, E (r), 



F> E (r) = A l (f ltE (r)+p l (k)g l , E (r)), 



r < r s , 
r > r s , 



(16a) 
(16b) 



where fi,E{r) and gi^ir) are respectively the solu- 
tions to the Schrodinger equation at energy E that 
are regular and irregular at the origin, and j3i(k) = 
— t&n(5i(k))/k 2l+1 is the scattering length coefficient. 
The irregular solutions are chosen in such a way, that 
as r — > 0, the ratio gi,E(r)/ fi,E(r) behaves as the ratio of 
the wave functions in free space ji(kr)/m(kr). The gen- 
eral form of the energy-dependent delta-shell pseudopo- 
tential follows as above, parameterized by some energy 



Vt, E Ar) = -f^^(3i(EQ)5(r-r s ) 



g'i,E ( r s) , d 



9i,E (r s ) dr 

(17) 

For the case of harmonic isotropic traps, with energy 
eigenvalues Eq = 2vq + I + 3/2, the Kummer function 
radial cigenfunctions Fi, E(l (r) = r l e~ r wi Ma (r) defined 
by Eq. (7), can be put into the canonical form, Eq. (16), 
by choosing 



/i,EoW =r l e 
9i,E (r) 

M{-v 



r M(-z/ ,/ + 3/2,r 2 ), 

„ + „ 2 [(2* + l)!!] 2 
21 + 1 
I -1/2,-1 + 1/2, r 2 ). 



(18a) 
(18b) 



Note that the pseudo potential depends on the trap- 
ping potential through the energy of the wave functions 
Eq = 2vq + 1 + ?>/2. Thus, the wave functions at any other 
energies but the bound state energies in the trap will ex- 
ponentially increase in the classically forbidden region. 
Since we evaluate the functions only at the shell radius 
r s in the construction of the pseudopotential, this is not 
a problem. Moreover, our pseudopotential, by construc- 
tion, only yields the correct asymptotic wave function 
at the energy Eq which do not have to be equal to the 
eigenvalue in the trap. To deal with this, we must set Eq 
equal to E self-consistently, as we discuss in the following 
section. 

To calculate the biorthonormal set of wave func- 
tions, we need the adjoint of the potential. Employing 
($1^4*) = (A T $|*), the adjoint potential is calculated to 
be 



V?{z) = a tan (Si) 



c 2 i5(r 



0) 



(19) 
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tions they are given by 

l/2A(l + A(%i,E(r s )///,i?(r- 5 )) 



1/2 + ci tan (<5 ; ) 



P> (r)=A i (/ i , E (r)+A(fc)^ E (r)). 



fi,E(r), 

(24a) 
(24b) 



FIG. 1: (a): The first (solid), second (dashed) and fifth (dash- 
dotted) wave function of the potential V = u5(r — r s )d/dr in a 
one dimensional infinite box of arbitrary length L, with r s — 
0.1L; (b) : the discontinuous wave functions of the adjoint 
potential, = d~ / dr(u8{r — r s ) .) . The wave functions form 
a complete biorthonormal set. 



This can be obtained easily from the eigenfunctions of 
the potential, Eqs. (16). 



IV. SELF-CONSISTENT SOLUTION FOR THE 
QUASI-HERMITIAN PSEUDOPOTENTIAL 



The boundary conditions for the adjoint set of wave func- 
tions, denoted P(r), are obtained in a manner similar Eq. 
(14), yielding 
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dr 



c\C2 tan 8i(k)P' 



= 0. 



(20) 

Because the adjoint potential contains a derivative of a 
delta function, the adjoint wave functions are discontin- 
uous at r = r s , as can be seen by integrating once more, 



1 



1 



-P < 



ci tan<5 i (fc)P < 



(21) 



For r > r s , however, these wave functions have the exact 
same behavior as F(r). 

As an example, consider two particles in an infinite 
spherical box of radius R, interacting through the po- 
tential V = u8(r — r s )d/dr. Integrating the Schrodinger 
equation over the potential yields the following boundary 
condition for the derivative of the wavefunction. 



dr 



9 T- 

or 



u — * 

or 



= 0, 



(22) 



The adjoint potential is given by d~ / dr(u8{r — r s ) ), 
where d~/dr signifies the derivative acting to the left. 
This leads to a discontinuous adjoint wave function, with 



0. 



(23) 



Both the wave function of the original potential and the 
derivative of the adjoint wavefunction are continuous, as 
can be seen in Fig. 1. Outside the radius of the shell 
the wave functions are identical, while the inside wave 
functions differ by their normalisation. 

The biorthogonal eigenfunctions obtained from Eqs 
(20,21) can be specialized to the case of atoms in a har- 
monic trap. In terms of the regular and irregular func- 



Associated with the pseudopotential, Eq. (17), is a 
biorthogonal set of eigenfunctions parameterized by an 
energy Eq. We seek the eigenfunctions for energies E 
that are the self-consistent solutions for two atoms in- 
teracting in the trap. Naively, one may attempt to use 
the self-consistent eigenvalues found from the Busch solu- 
tions, Eq. (9), as the self-consistent parameter Eq. These 
resulting wave functions do not, however, give rise to a 
biorthogonal set since the regularization will vary for each 
eigenfunction. Instead, we must fix the parameter Eq in 
Eq. (17) at an arbitrary value to find a given biorthog- 
onal set, and then find the value of Eq that allows us to 
satisfy all the boundary conditions in the trap. A side 
effect of this procedure is that for energies E ^ Eq, the 
scattering phase shift arising from the pseudopotential 
is not equal to the scattering phase shift of the true po- 
tential at that energy. We are thus faced with the task 
of ensuring self-consistent scattering phase shifts at the 
same time as we satisfy the trap's boundary condition. 
We achieve this through a dual parameterization as dis- 
cussed below. 

Consider the solutions at an energy E of the pseudopo- 
tential delta-shell parameterized by an unknown energy 
Eq. Given the solution outside the shell in canonical 
form, Ai(fi(E) — f3gi(E)), and the boundary conditions, 
Eq. (14), the scattering length function is, 



ME,E ) 



(fi9io-fl9io)Pi(Eo) 



(fto/Mfi g[ - fl 9i) + (9w g[ - gi 9'MEq) ' 



(25) 



where the regular and irregular functions, /; and gi are 
evaluated at the shell radius, r = r s . Here, giQ — gi t E i 
9i = 9i.E, fis) = fi,E , fi = fi,E, and g' m = g' l Eg denotes 
the radial derivative. The scattering phase shift that is 
imposed on the wave function by the pseudopotential is 
strongly energy dependent, matching the true interaction 
potential only when E = Eq, so that (3i(E = Eq,Eq) = 
Pi(Eq). To better understand the behavior of the scat- 
tering length coefficient, consider the case of free-space 
collisions where Ji e = ji{kr s ) and g^E = ni{kr s ). In the 
limit as r s — > 0, the s-wave and p-wave scattering length 
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FIG. 2: The dressed scattering length (3(E,Eo) that is im- 
posed on the wave function by the pseudopotential as a func- 
tion of the energy, for different bare scattering lengths /3(E) 
and E = 1. Solid line: 13(E) = -5, dashed: /3(E) = 1; 
dashed-dotted: /3(E) = 10. For E = Eo the bare and dressed 
scattering lengths agree, /3(E,E ) — /3(E), and for other en- 
ergies there is no agreement. For fixed Eo the position of the 
pole depends on the bare scattering length of the interaction 
potential, while the width of the pole depends on the shell ra- 
dius r s . All units are normalized to the characteristic scales 
of the harmonic oscillator as discussed in the text. 



coefficients that follow from Eq. (25) are, 
k/3 (E ) 



Po(E,E ) 



0i(E,Eo) 



k + r s (kl-k 2 )fo(E y 

3fc 3 /3!(^ ) 

2kl + kkl-f- (fc 2 -P)/M£o)' 



(26a) 
(26b) 



Note the existence of a pole in (3 at an energy depending 
on the shell radius. In the case of s-waves, the position 
of the pole goes to k = oo as r s — > and we recover 
/3o{E,E ) = (k/k )P {Eo) for finite energies. Thus, for 
s-waves, one can construct a contact potential simply 
using the bare energy-dependent scattering length. In 
contrast, for p-waves, the limit r s — > is singular. The 
pole occurs at k = ko and the scattering length is inde- 
terminate. One requires, therefore, a pseudopotential for 
higher partial waves with finite radius r s . For the case of 
two atoms in a harmonic trap interacting via p-wave col- 
lisions, the outside wave function is given by the Kummer 
U function. Eq. (25) then reduces for small r s to, 



Pi(E,E ) 



3(3i(E ) 



E)P!(E ) 



(27) 



The width of the pole depends inversely on the shell ra- 
dius r s . An example of /3i(E,E ) with fixed E can be 
seen in Fig. 2. 

We can now calculate a self-consistent solution for in- 
teracting atoms in the same trap. For a given interaction 
potential we calculate (analytically where possible, nu- 
merically otherwise) the bare scattering length functions 
Pl(E ) = -tan(8[(E ))/kl l+1 . Using Eq. (25) we find 



the dressed scattering length function of the pseudopo- 
tential (3i(E, Eo). The generalized Busch eigenvalues are 
then found from Eq. (9), with /3 t (E) -> /3i(E,E ), 



Pt(E,E Q ) = (- 



TT 



r(Z + 3/2) 



(2Z + 1)!! 



T(-v -1-1/2) 



(28) 

The solutions of this transcendental equation with E = 
2v + I + 3/2 can be found numerically for arbitrary 
Eo- Given the energy eigenvalues, one can determine 
the biorthonormal basis, parameterized by E , using 
the Kummer equations and boundary conditions Eqs. 
(14,21). For two atoms in the same trap, the self- 
consistent solution follows by setting E equal to these 
eigenvalues. In the next section, we solve the problem 
of two interacting atoms in separated traps via this self- 
consistent procedure. 



V. TRAP-INDUCED RESONANCES FOR 
HIGHER PARTIAL WAVE SCATTERING 

We consider the application of our formalism to study 
the interaction of two atoms, trapped separately in dis- 
placed harmonic wells. We assume each well is seen 
only by one of the distinguishable atoms, as in the state- 
dependent traps associated with a polarization gradient 
lattice [1]. The Hamiltonian of the system is 



1 



2 2 

TT Pi , P2 1 ± 2 

H = h h —TflLO 

2m 2m 2 



ri + 



Az 



+ 



-muj 



Az 



r2 - 



^nt(ri-r 2 ), (29) 



where Az is the vector separating the traps and V- ln t (ri — 
Y2) is the true interaction potential between the atoms, 
which we will model using our pseudopotential. Due 
to the quadratic form of the Hamiltonian, one can sep- 
arate the noninteracting center-of-mass from the rela- 
tive coordinate. With our choice of units, i.e. h = 1, 
H = m/2 = 1, w = 1 and the characteristic length 
zo = y^h/JJiLuj, this is governed by the following Hamil- 
tonian. 



Ho 



H iel = H + H(Az), (30a) 

1 d 2 1(1 + 1) 2 \ „ , , ,„„, . 

■-^ + A 72-^ + ^ 2 j+^tW, (30b) 

H(Az) 



Azrcos6»+ -Az 2 , (30c) 



where 9 is the polar angle between the relative coordinate 
and Az. We have studied this problem previously in 
the case of s-wave collisions and observed the existence 
of trap-induced resonances when the potential energy of 
the trap shifts a molecular bound state of the interaction 
into resonance with a vibrational state of the trap [11]. 
Here, we study the generalization to higher partial wave 
scattering. 
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FIG. 3: Partial wave scattering lengths for the square well 
potential as a function of energy in harmonic oscillator units. 
Dashed line: s-wave scattering length (in units of the char- 
acteristic length of a harmonic oscillator, zo); Solid line: the 
scattering volume (3(E) for p- waves in units of Zq. 



A numerical solution follows from the quasi-Hcrmitian 
pseudopotcntial construction. We can construct a com- 
plete biorthogonal set from the eigenfunctions of Hq with 
Vint set to the delta-shell pseudopotential, parameterized 
by an unknown energy E and scattering phase shift 
Si(Eo) according to the procedure in Sec. (IV). We use 
these functions to find a matrix representation of the per- 
turbation H(Az). Diagonalizing yields the total eigen- 
value E. The self-consistent solution is the one where 
E = E + Az 2 /2, so that E is the kinetic energy of the 
collision. 

We test this procedure with a toy model of the molec- 
ular binding potential: a spherically symmetric step po- 
tential, of depth Vq and radius Ro, whose scattering phase 
shift can be found analytically. A plot of the scatter- 
ing length can be seen in Fig. 3. We compare the self- 
consistent solutions based on the pseudopotcntial with 
those generated from an orthonormal basis - eigenstates 
of H with Vi n t taken to be the exact potential. The lat- 
ter basis is determined by matching Kummer function 
radial wave functions across the step radius as in [19]. In 
both cases, we neglect scattering beyond p-wave. We only 
show states with angular momentum projection m = 
along the direction of the trap separation. Due to the 
axial symmetry of the problem, the states with different 
to are not coupled to each other. 

The agreement is excellent as can be seen in Fig. 4. 
The plot shows some interesting features of trapped sep- 
arated atoms. For our step- well, there is a p-wave bound 
state at energy E = — 2; the s-wave bound state is very 
deep and not shown. The p-wave bound state is shifted 
up quadratically by the trap, leading to trap-induced 
resonances where it crosses the higher lying vibrational 
states, due to the mixing of the partial waves at finite 
separation described below. Note that it only mixes with 
one of the states of each degenerate subspace, leading to 
both an avoided and an actual crossing for each of the 
degenerate manifolds. At Az = 0, the higher lying states 
are p, d, s, f and p states, respectively, as noted in 




FIG. 4: Eigenspectrum for two atoms in displaced harmonic 
oscillators as a function of trap separation, for the example of 
a square- well interaction potential with a radius of Ro = 0.1 
and depth Vb = 489.9 in harmonic oscillator units. Only the 
states with angular momentum projection in the direction of 
the separation m = are plotted. Interaction between higher 
partial waves than p-waves were neglected. The slash-dotted 
line is the solution based on the biorthonormal eigenfunctions 
of the pseudopotential, while the thin line shows the numerical 
solution to the exact problem. Each state is denoted by its 
angular momentum I for no trap separation. For finite Az the 
states are mixed due to the separation. Several trap-induced 
resonances due to the shifted molecular p-wave boundstate 
are visible. The agreement is excellent. 



Fig. 4. Only s and p states are shifted in our approxi- 
mation. For finite separation, the vibrational states are 
superpositions of partial waves with respect to the in- 
teraction center. In each subspace of vibrational states 
that has a set of degenerate partial waves with respect to 
the trap center without the interaction, there is a linear 
transformation of the degenerate eigenspace such that 
the j th eigenvector contains no partial wave components 
for 1 = 0... j—2. Thus, as seen in our previous work, [11], 
with only s-wave interactions, the second excited vibra- 
tional manifold contained one state which is unshiftcd. 
In our case, this state is affected by p-wave scattering at 
finite separation. In general, including up to k partial 
waves in the interaction, k states within a degenerate 
subspace will be shifted. A small avoided crossing due 
to higher-order interactions can be seen between these 
states at a separation of 0.9. Similar behavior is seen for 
the third vibrational subspace, with an avoided crossing 
between the sixth and seventh state at a separation of 
Az = 1.3. Mixing of partial wave solution is generally 
seen in anisotropic traps [16]. 
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VI. SUMMARY 

In this paper, we considered a quasi-Hermitian pseu- 
dopotcntal for higher partial wave scattering based on 
a delta-shell potential with energy-dependent regulariza- 
tion. In contrast to all other generalizations for I > 0, 
the regularization derived here contains only a first order 
derivative, rendering it quasi-Hermitian and supporting 
a complete biorthogonal set of eigenfunctions. Such a 
set is useful for studying interactions of atoms in tightly 
confining traps as it forms a basis for easily diagonaliz- 
ing the relevant Hamiltonian. Because the regularization 
is energy dependent, we must solve such problems self- 
consistently, so that the eigenvalues of the diagonaliza- 
tion match those of the parameters in the pseudopoten- 
tial. 

We have tested our construction by studying the prob- 
lem of two atoms trapped in separated wells, with a 
square well to model the inter-atomic interaction. The 
solutions obtained from our pseudopotential show excel- 
lent agreement with the numerical solutions to the com- 
plete problem. We also showed that a p-wave bound state 
close to dissociation leads to trap-induced resonances, 
which have been calculated earlier for s-wave scatter- 
ing. This is particularly interesting for identical fermions 
which cannot interact via s-wave collisions due to their 
quantum statistics. 

In future research, we plan to apply our method to 
study realistic atomic systems. A particular relevant 
case is cesium collisions. Due to a bound state of the 
cesium dimer almost exactly at dissociation, the collision 
properties are anomalous. We have already seen large 
trap-induced resonances for s-wave scattering in previ- 
ous theoretical studies [12]. Because cesium has large 
spin-relaxation, high partial wave scattering can play an 
important role. Trap-induced resonances should then be 
exhibited in multichannel scattering processes. 



the kinetic energy operator when its domain is extended 
to include the space of the irregular functions at the ori- 
gin [31]. Surprisingly the non-Hermitian parts of these 
operators cancel exactly, as can be seen as follows. The 
antihermitian part of our Hamiltonian is 



f 

Jo 

1 

2 



(F*HF j 
d 



F HF?) dr 
d 



(rFnlirF^+ajFr^rFj) 



(A.l) 



Here. 



and Oj arc the energy-dependent scattering 
lengths at the eigenenergies 2z/j + 3/2 and 2vj +3/2 
for Fi and Fj, respectively. The wave functions are 

Kummcr {/-functions, given by Ce~ r (M(— Vi, 3/2, r 2 ) — 
a,l/rM(— Vi — 1/2, 1/2, r 2 )), where C is a constant. 

2 

Pulling out the factor of e~ r and using the fact that 
M(a, b, r) = 1 + ar/b + wc obtain 



HAij - aj +di- a t + 



+ a,j - a,i - 



= 0. (A.2) 



Since this is zero for all the Hamiltonian for the s- 
wave pseudopotential is Hcrmitian. Thus the eigenfunc- 
tions form an orthonormal set. No such cancellation can 
occur for higher partial waves, as discussed in Sec. Ill 
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